Why does the Unruh effect rely on Lorentz invariance, 
while Hawking radiation does not ? 
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We show that without Lorentz invariance, the Unruh effect does not exist. We use modified 
dispersion relations and describe in turn: the non-thermal nature of the vacuum (defined in the 
preferred frame) restricted to the Rindler wedge, the loss of the KMS property of the Wigthman 
function, the transition amplitudes and transition rates of a uniformaly accelerated detector. This 
situation seems to contrast with the Hawking radiation of acoustic black holes, which under certain 
assumptions has been shown to be robust to a breaking of Lorentz symmetry. We explain this 
discrepancy. 



It is well understood that Unruh effect and Hawking 
radiation are physically distinct phenomena but with a 
common root, if not a physical at least a mathematical 
one [1]. Yet, the reason of this analogy is not fully un- 
derstood. 

The situation seems in fact paradoxical. On one hand, 
we are going to show that the Unruh effect is inseparable 
from Lorentz symmetry: without the latter, the former 
does not exist. On the other hand, the study of black 
holes analogues in supersonic fluids and Bose-Einstein 
condensates, the so-called acoustic black holes or dumb 
holes, has shown (under certain assumptions) that the 
defining properties of Hawking radiation (in that case 
a stationary and thermal flux of phonons escaping from 
the sonic horizon to infinity) are robust to a breaking of 
Lorentz symmetry, sec [2] for a review. So how can they 
be related ? A reconciliation of these two results seems 
necessary in view of the analogy mentioned above. 

The fundamental role played by Lorentz symmetry in 
the existence of the Unruh effect is particularly clear 
from the algebraic proof of this one [3, 4]. This proof 
establishes, at the level of the operator algebra, that the 
Minkowski vacuum restricted to a Rindler wedge is a 
thermal state with respect to the boost parameter. This 
theorem is proved in the framework of axiomatic field the- 
ories, which makes the instrumental role of Lorentz sym- 
metry abundantly clear. Yet one can wonder whether, in 
the absence of the invariance under boosts, the Unruh ef- 
fect exists in a more operational sense (outside transients 
of course). As we will show, the answer is negative. 

We conducted a complete analysis of the Unruh ef- 
fect without Lorentz invariance. We examine the map- 
ping between Minkowski and Rindler Fock spaces (Unruh 



modes and Bogoliubov transformations), the Wightman 
function, and the transition amplitudes and transition 
rates of a detector coupled to the field. Lorentz symme- 
try considerably constrains the structure of field theories 
and renders all these descriptions of the effect equivalent. 
The unifying role of Lorentz symmetry and the instru- 
mental role of the stable vacuum are recalled in sec. I. 
Without Lorentz symmetry, each quantity answers to a 
different question as will be explained in sec. II. Finally 
we will compare Unruh and Hawking effects in section 
III. The reader will find the main results summarized 
and discussed in the subsection E of each section. 



I. UNRUH EFFECT WITH LORENTZ 
INVARIANCE 

We set up our notations and review the Lorentz invari- 
ant case. The contend of this section is not new but our 
presentation, which lays emphasis on the role of Lorentz 
symmetry, may be original, see in particular sec. IE. We 
refer to the review articles [5, 6] for more details. 

We consider the simple case of a free massless scalar 
field f(t,x) in Minkowski space-time. In incrtial frames 
equiped with a global and cartesian coordinate system, 
i.e. ds 2 = —dt 2 + 5ijdx l dxi with units c = 1, the wave 
equation is 

(d 2 -^0^)^ = 0. (1) 

The field is assumed to be neutral so its creation and 
annihilation parts are conjugate from one another 

¥>(i,x)= d 3 k (a k </? k (i,x) + 4<^(t,x)) , (2) 



where 
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are the plane wave solutions of (1) with positive fre- 
quency w.r.t the inertial time coordinates, i.e. idt^k = 
oJkfk > w fc > 0. For masslcss fields the dispersion relation 
is of course uj k = k = |k|. The state of the field is the 
Lorentz invariant ground state defined by <Xk|0M) =0 for 
all k. 

A two-level point detector is coupled to the field as 
described by the interaction Hamiltonian 

H int =g(e iEr \+)(-\+e- iET \-)(+\) <p[x a (r)]. (4) 

g is a dimensionless coupling constant. are respec- 
tively the ground and excited energy eigenstates. They 
are separated by the energy gap E > in the detector's 
rest frame, r is the detector's proper time and x a (r) the 
detector's trajectory. The latter will be constrained to 
be linear, and we can choose the coordinates such that 
the detector moves along the z-axis. Inertial (In) timelike 
trajectories with velocity /3 are therefore parametrized by 

Inertial : [ 7 r, 0, 0, 7/3r] , (5) 

where 7 = (1 — /3 2 )~ 1//2 is the Lorentz factor. Uni- 
formly Accelerated (UA) detectors with acceleration a — 

/ 2 a 2 b \ 1 

yiab ^dpr^dpr) follow the trajectory 

UA : [a -1 sinh(ar), 0, 0, ± a -1 cosh(ar)] . (6) 
The corresponding Lorentz factor is 

T(r) = ^ = cosh(ar) . (7) 

The literature on the Unruh effect is sometimes con- 
fused on the definition of the effect. There are two dif- 
ferent aspects to it. The first aspect is the thermal 
nature of the Minkowski vacuum for observables with 
support in a Rindler wedge. The R(ight) and L(eft) 
Rindler wedges of Minkowski space are the open sets 
{(t,x,y, z)/z ^ 0, |t| < ±z}. They are both static, glob- 
ally hyperbolic space-times, so that a consistent quantum 
theory can be defined on them [7]. This theory is not 
equivalent to the usual quantization in the full Minkowski 
space-time. Rather, the Minkowski vacuum corresponds 
to a thermal state of Rindler quanta at "temperature 
l/27r" (in arbitrary units). The second aspect is the op- 
erational meaning of the previous result in terms of the 
response of a UA detector [8]. (As we can see from (6), 
a UA trajectory is confined to a Rindler wedge.) To wit, 
the S'-matrix elements and transition rates of the detec- 
tor verify detailed balance. 

The thermal nature of the Minkowski vacuum re- 
stricted to a Rindler wedge is recalled in section I A and 
the Wightman function is described in IB. We then con- 
sider the dynamics of the detector introduced at eq. (4). 
We will see respectively in sec. I C and I D that the Bo- 
goliubov coefficients of section I A are proportional to 
the S'-matrix elements of the processes |=f) — > |±) + k, 
and that the Wightman function transmits its properties 



to the transition rates which are essentially the Fourier 
transform of the former. We will insist on how the re- 
sults are related by Lorentz symmetry. We therefore do 
not choose to define the Unruh effect by one or the other 
aspect because by Lorentz invariance they are equivalent. 

A. The Minkowski state for observables in a 
Rindler wedge 

The proof begins with the definition of the field theory 
in the static, globally hyperbolic Rindler wedges. The 
metric in the R wedge can be brought in the static form 
by the following change of coordinates 

t = K~ 1 e K< ' sinh(rer7) , z = K~ 1 e K( ' cosh(ray) (8) 

where re is an arbitrary energy scale. We choose units re = 
1. The line element is ds 2 — e 2C -(—dr] 2 +dC, 2 ) + dx 2 + dy 2 . 
The timclike Killing vector corresponds in the R wedge 
to the generator zd t + td z of boosts in the z-direction. 
The future and past horizons are respectively 

^ ± = {V -> ±00, C -> -00, 77 ± C finite} . (9) 

The null coordinates u = r\ — ( and v = 77 + £ will be 
usefull. Letting xj^ = (x,y), the field equation (1) in 
these coordinates is 

(d 2 d 2 e 2 < (d 2 x + d 2 )) pfa, C, x ± ) = . (10) 

Since the metric is static, the solutions of (10) can be 
classified according to the eigenfunction of the timelike 
Killing vector field id n — A where A > 0, given by 



Similarly in the L quadrant we can introduce the coordi- 
nates (77, x, y, C) defined by 

t = — sinh(r7) , z = — cosh(77) (12) 

We chose dt/dfj < because the boost Killing vector field 
td z +zdt is timelike and past-directed in that wedge. This 
convention implies that the modes defined by 

^(r7,C,xx) = ^JrU,xx), A>0. (13) 

are positive frequency Rindler modes. This double 
familly of modes form a complete orthonormal basis on 
which the field can be decomposed 

¥>(t,C,xl) = J d\ Jd 2 k ± (a^ k± <p^ k± +at k± vti ± +h.c.) 

(14) 

We designate by 10^) the vacua in the R and L wedges 
respectively. 
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To an observer living in the i?-wedge, that is an ob- 
server measuring observables Or <£> 1l, the Minkowski 
vacuum appears to be the mixed state p defined by 

Tr = Tr (\0 M )(0 M \6 R ® 1l) (15) 

for every observable, that is 

P = Tr i (|0 M )(0M|) (16) 

where the partial trace Tr^ is over the Hilbert space of 
the theory defined in the L-wedge. If the field theory is 
free, it is possible to calculate p explicitcly. To this end, 
one first establishes the unitary map between |0m) and 
the Rindler Fock states, or equivalently, since the theory 
is free, a unitary transformation between the positive fre- 
quency Minkowski modes (3) and the Rindler modes (11) 
and (13). From now on we do not write the subscript kj^ 
anymore since it is obvious that the map does not mix 
different transverse wavevectors. 

This map is found by appealing for the stability of 
the Minkowski vacuum and the concomitant analytic 
properties of the Minkowski modes. Namely, the sta- 
bility of the vacuum in any inertial frame is tanta- 
mount to the analyticity and boundedness of the posi- 
tive frequency Minkowski modes (3) in the domain T = 
{x + iy\y° < 0,\y°\ > |y|}- Since any linear combination 
of positive frequency Minkowski modes is also analytic 
and bounded in T, the idea is to define a new familly 
of modes <Pq, defined in the complete Minkowski space, 
such that i) they are eigenmodes of the boost Killing 
vector zd t +td z , hence of the form (no Rindler-frequency 
mixing) 

fl>0, v>K = an^ + Ai¥»n* (17a) 
fl<0, <p% = a m <pf n] +P ln \<P\h*\ , (17b) 

ii) admit the decomposition into positive frequency 
Minkowski modes 

/+oo 
dk z 4 A p£ , (18a) 
-co 



/+oo 
-co 



dk z B nMz ip k 



M * 



(18b) 



iii) and are orthonormal and complete (w.r.t. the Klein- 
Gordon product). This basis of solutions of the wave 
equation, called the Unruh modes, characterizes the 
Minkowski vacuum and contains within its definition (17) 
the mapping we are looking for [8] . 

For the function ip^ introduced in eqs. (17) to be de- 
fined on the entire Minkowski space, we need to extend 
the definition of the Rindler modes ip R to the L-wedge. 
From now on ip R denotes the function equal to the r.h.s. 
of (11) in the R- wedge, and equal to zero in the L-wedge. 
We define similarly <p\ the extension of the L-Rindlcr 
modes (13) to the i?-wedge. 

The task to find the Bogoliubov coefficients an and /3q 
such that ifQ as defined by eqs. (17) verifies condition ii) 



is simplified by the following theorem [1]: any solution 
of the Klein-Gordon equation in dimensions more than 
2 is characterized by its restriction to either H + or H~ . 
Written in terms of the null coordinates U = t — z and 
V = t + z, the positive frequency Minkowski modes 



M 

Vk z « e 



i(u k -k z )V/2-i(u k +k z )U/2 



(19) 



are analytic and bounded functions for the complex val- 
ues of U and V such that Im(V) < and Im(f7) < 
because uj\ = k 2 z + k\ > k z . We therefore require that 
the restriction of the linear combination (17) on, say the 
future horizon shares the same property. The asymptotic 
form of (11) near % + is the sum of two terms 



(fc±/2) tA 

r(i + *A)' 



(fc±/2)~ iA 
r(l-iA) 



(20) 



We do not need the expression of N\ for the moment. On 



H- 



oo and v is finite, so the first term is singular. 



When the modes are superposed to form wave-packets, 
the term which rapidly oscillates as the wave-packet nears 
the horizon does not contribute. The reason for consid- 
ering wavepackets is that fields are operator valued dis- 
tributions, and therefore must be smeared. (We did not 
bring up this issue before because it is the only place in 
this section where it is mandatory to use wavepackets 
instead of modes.) So we retain only the second term. 
Since the null coordinates are related by V\r — e v and 
V\ L = -e c , the r.h.s. of (17) is 

O>0,^ ~ an\V\- in (e(V) + 6(-V)^) ,(21) 

and similarly for Q < 0. The unique analytic continua- 
tion of V~ in which is bounded in Im(V) < is (V— ie)~' lU 
where the branch cut of the logarithm extended to the 
complex plane is chosen along the negative real axis. The 
relative weight is therefore 



£J = e -7r|n| 
an 



and (Pq is normalised if 



an 



1 

I _ e -27r|0| ' 



(22) 



(23) 



an can be chosen real. 

A direct calculation finally shows that the (p^ form a 
complete familly, so that the field can be represented as 
follows 



= TdSl fd 2 k ± (a^ + h.c.) (24) 

J — oo J 



and the Minkowski vacuum is characterized by 

ag|0 M )=0. (25) 
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One can invert the relations (17) to relate the creation 
and annihilation operators 

Q>0, ag = a n ag- faatf ^Ua^U' 1 , 

n<0, ag = ainiafo - 0\ci\a$ = Uaf^U' 1 , (26) 

where U is the squeezing operator 



u = n° x p { ar § th (^) («S a N 



Lt -4i4iU ( 27 ) 



This is the unitary map we where after. The Minkowski 
vacuum is indeed related to the Rindlcr Fock states by 

|0 M ) = U\0 R )\0 L ) 

oo 

= n E e " n27r|n| i n «)i^)- ( 2g ) 

f!,kj_ n=0 

In consequence, an observer living in say the R wedge 
interpretes the Minkowski vacuum as a thermal bath of 
i?-quanta at the temperature T — l/2n 



P= II e- n ^ n \\n R ){n R \ 



(29) 



n,kj. 



We finally mention that the Bogoliubov coefficients 
have two different physical interpretations. The first is 
provided by (26) : they describe the mean values and cor- 
relations of Rindlcr quanta in the Minkowski vacuum 



(0 M \a^a R ,\0 M ) = 8(\- \')\/3 x \ 2 , 
(0mK4|0 m > = S(X-X')a x p x . 



(30) 



The second interpretation is given by their proportional- 
ity to the S'-matrix elements, see eq. 50. These identities 
as well as the expressions of the coefficients A and B in 
eq. (18) are established in appendix A. 



B. Wight man function 

The interpretation of \0m) either as the ground state or 
a thermal state is reflected in the corresponding expres- 
sions of the two-point Wightman function. For arbitrary 
events with inertial coordinates y and x + y, 



W(x) 



{0 M \tp{x + y)<p(y)\0M) 

- 1 . J dk g-i^x ^,ik|x| _ e -*fc|x|^j e -ek 



8tt 2 |x| 

1 



1 



47T 2 (x° - it) 2 - x 2 ' 



(31) 



The last expression exhibits the analyticity of W in the 
lower half complex x°-plane and can be obtained from 
the integral representation by the introduction of a regu- 
lator e~ cfc . This regulator implements the stability of the 
vacuum since, as we already saw, the Minkowski modes 



(3) are analytic and bounded in the complex space-time 
domain T = {x + iy\ y° < 0, \y°\ > |y|}. By Lorentz in- 
variance, this in turn is equivalent to say that the one 
dimensional section of the Wightman function on the in- 
ertial straight lines is an analytic function of the com- 
plexified affine parameter t of the geodesies in the lower 
half plane Im(r) < 0. The expression of the Wightman 
function, when both points belong to a common inertial 
trajectory (5), 



W/„(T) = (0 M Hxi n (T + T')]ip[x In (T')}\0 M ) 
1 1 



47T 2 (r - ie) 2 



(32) 



is indeed an analytic function in the lower half complex 
r-plane, where r 2 = (x ) 2 — |x| 2 . 

The expression of the Wightman function when x and 
y are arbitrary points in the R wedge is not illuminat- 
ing. It becomes interesting only if the points are on an 
orbit of the generators of boosts such as the linearly UA 
trajectory (6) 

W ua (T) = (0M\(p[Xua{T + T')]ip[x ua {T')}\0 M ) 

= _^ 1 (33 ) 

16tt 2 sinh 2 (f -ie) 1 ' 

This function enjoys two important properties. First, it 
depends only on the difference r of the proper times, 
because the vacuum is Lorentz invariant and a shift of t 
corresponds to a boost. Second, (33) verifies the KMS 
condition 



t + i- 



2n 



W ua (-T) 



(34) 



which is the definition of an equilibrium state of temper- 
ature a/2-K. 

Finally, the expression of the Wightman function on 
two points in opposite Rindlcr wedges, both at proper 
distance a -1 from the horizons, that is the points 



t R = — sinh ar , 
a 

zr = - cosh ar , 
a 

xij; = 0, 



ZL 



sinh ar' . 



cosh ar' . 



x ±L =0 



(35) 



is given by 



W RL (6) = ^mHxr{t)]^[x l {t')]\Q m ) 

„2 



167T 2 cosh (aS) 



(36) 



where 6 = (t— t')/2 (recall that — t' is the future directed 
proper time on the UA trajectory in the L- wedge), see 
also sec. IV. C in [9]. The correlations are maximal for 
pairs of conjugate points r = r' and decay exponentialy 
outside the region |<5| > I /a. The Fourier transform of 
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the Wightman function w.r.t. S is 

/+oo 
dSe iX5 W RL (S) 
-oo 

_JA| e^ x l 2a 
~ 8tt 1 - e ~ vX / a ' 



One recognizes the factor ot\f3\ of the Bogoliubov coef- 
ficients (22) and (23). To show this directly, we use the 
expansion of the field in terms of Unruh modes (24) and 
evaluate the expression on the points (35) with the help 
of the relations (17) 

/r+oo 
d 2 k ± j d\a x (3x<Px{x R )<p^{x L ) 

d?k±J d\e- taX( - T - T > \K lX (k ± )\ 2 . (38) 



amplitudes are literaly the expression of the stability of 
the vacuum. 

For the UA trajectory, we choose the scale k = a in 
the coordinate system (8) so that the time coordinate r\ 
coincides with the proper time along the UA trajectory 
(37) (6), which now coincides with the coordinate curve £ = 0. 
We introduce the shorthand notations 

k± = oj k ± k z (44) 

We look for an analytic expression of the integral in 

-ig/a 



A 



UA,LI 
±,k 



/+oo 
dx e ±i%*+i(k-e°-k + e-*)/2 (45) 
-oo 



The identity 



ttA 



POO 

/ duu\K iX {u)\ 2 = — - 
J Q 2smh(7rA) 



(39) 



combined with the expressions (23) and (22) finally gives 



Wrl(S) 



+ 00 



XdX e -^ x ^ a 



47T 2 1 - e -27r|A|/o 



-ia\X\(r-T') 



(40) 



whose Fourier transform w.r.t. to S = (r — t')/2 is (37). 



C. S-matrix elements 

The S-matrix element of the process |=f) — > |±) +k, i.e. 
the exitation (desexitation) of the detector accompanied 
by the emission of a particle of momentum k from the 
Minkowski vacuum |0m) are at the lowest order 

A L J M = (OmK ® (±| f e ~ 4 |t) ® |o M ) 



-w 



^2^(277)3 

For the inertial trajectories (5) we get 



/+oo 
^ T e i[±BT+u t f(r)-kx(r)] _ 
-oo 



where 



t In, LI 
*±,k 



-«5 



*(u/ fc ±25) 



7 (u; fe -/3fc z ) > 



(42) 



(43) 



is the energy of the scalar quantum in the rest frame of 
the detector. This quantity is stricktly positive in any 
Lorcntz frame, as a consequence of the stability of the 
Lorentz invariant vacuum. Hence the amplitude to spon- 
taneously emit a quantum from the ground state van- 
ishes. The Dirac distribution S(ui' k — E) in the amplitude 
A- is the expression of the conservation of energy. These 



I = k 2 z + > fc 2 , so that both k± are 
positive. We assume first kj ^ ^ 0. After the change of 
variables y = x + In sjk-jk+ the amplitude reads 



A UA ' LI oc 



±iE/2a 



+oo 



dye 



±i^- ikj_ sinh(y) 



(46) 



We evaluate it as a contour integral along the rect- 
angle with edges on Im(y) = and Im(y) = — it/2. 
The integrals along the vertical axis vanish and the 
integral along Im(y) = — 7r/2 is the integral repre- 
sentation of the modified Besscl function 2K u (z) = 



A 



UA,LI 
±,k 



-i2g 



v/2^(2^) 3 a 



K. 



a 



• (47) 



The case kj^ = must be treated separately but the ex- 
pressions turn out to be the limiting values of the r.h.s. 
of (47) [10]. Since K_ v (x) = K v (x), the corresponding 
probabilities P±m = \A l ^ LI \ 2 differ by a Boltzman fac- 
tor 



+ ,k _ e -2nE/a 



(48) 



The origin of the Boltzmann ratio and the relationship 
with the analytic structure of the Wightman function are 
perhaps better understood from the saddle point approx- 
imation of (46). The saddle points of A± are 



s+ = argch(S//cj 



in2ir . 



x n = iir + x+ . (49) 



First, the periodicity of these saddle points reflects the 
periodicity of the poles of the Wightman function. (This 
should not be a surprise because the Wightman function 
and the amplitudes are integrals of the mode functions) . 
Second, once exponentiated, the relative shift in between 
the positions of the saddle points gives the Boltzmann 
factor e - iEx ~/ a = e *E/a e -iExi/a^ ( A more comp i ete 
analysis can be found in appendix C and section II D). 

Last but not least, the S- matrix elements are propor- 
tional to the Bogoliubov coefficients between Minkowski 



6 



and Rindler modes (these expressions are shown in ap- 
pendix A), 



A 



UA,LI 



■ 9 ts f k ±- 

-l-K-i\/a 

7r 1 \ a 



sinh(7rA/a) 



a 



a ua,li_.9 k f k ±\J sinh(7rA/a) 

A +M - *~ K iX/a — ) Y P\,k z M± ■ ( 50 ) 

This is the consequence of the facts that trajectories of 
UA observers are curves of constant Rindler coordinate 
C, and along these trajectories r = ai] is proportional 
to the Rindler time coordinate. In consequence the am- 
plitudes and Bogoliubov coefficients are essentially given 
by the same integral. Lorentz symmetry thus endows the 
Bogoliubov coefficents with a dynamical interpretation. 



D. Transition rates 

From the S'-matrix elements of the previous section one 
forms the inclusive probabilities 



P± = d 3 k \A 



±,k| 



(51) 



by summing over the final states of the field. They can be 
calculated directly from the expressions of the probability 
amplitudes given in the previous section, or alternately 
by exchanging the order of integration over k and time. 
The calculation is facilitated by keeping the upper bound 
of the time integration in (41) finite and taking the limit 
r — > oo at the end. Calculations with the first method 
can be found in [6]. We adopt the second method [5] 
which relates the probability to the Wightman function 
of the field, 

/T r+co 
dn / drae^^Wfa.-n-Ta) (52) 
-oo JO 

where the Wightman function is evaluated at two points 
on the trajectory of the detector. The corresponding 
transition rates of the detector are defined by 



R±(r) = 



dP ± 

~dV 



r-t-oo 

= 2g 2 Re / dr' e TtEr ' W(t, t-t'). (53) 
J o 

From the expressions (32) and (33), we see that 
W(t,t — t') in the integrand depends only on the dif- 
ference of its arguments r' both for inertial (5) and UA 
trajectories (6). The rates R± are thus time-independent. 



Again, this is because both trajectories are orbits of 
a Killing vector of Minkowski space and because the 
Minkowski vacuum state is annihilated by the corre- 
sponding generators. Namely, inertial trajectories are 
invariant by time-translation since d T \i n — jd t , and for 
UA trajectories a translation along the proper time is a 
boost since d T \ ua oc td z + zd t . 

With this simplification, one can write (53) as an inte- 
gral along the entire real line 



i?4 J = 9 Z Re 



/+oo 
-OO 



W(t) 



(54) 



where W(r) is either (32) or (33). One can calculate this 
integral by the method of residues. On inertial trajec- 
tories, the Wightman function (32) has a double pole a 
t = ie and is analytic in the lower half complex r-plane 
(since e > 0), which yields 



R T + n ' LI = , 



R 



LI,In 



2tt 



(55) 



The first rate vanishes because of the stability of the vac- 
uum in all Lorentz frames. 

For UA trajectories, the Wightman function (33) has 
now a countable family of double poles t„ = i(n2ir/a + e) 
periodically spaced on the imaginary axis. As a result 



R 



UA.LI 



= ± 



2tt 



1 



1 



(56) 



As for the amplitudes, the Boltzmann factor follows di- 
rectly from this periodicity and the stability of the vac- 
uum (the ie prescription says that the pole To counts in 
R- but not in R + ) 



R 



UA.LI 

+ 



R 



UA,LI 



-2ttE/ci 



(57) 



E. The double role of Lorentz symmetry 

Lorentz invariance is instrumental in the previous re- 
sults in two respects. First, to ensure the stability of 
the ground state in every frame. The latter is sole re- 
sponsible of the properties of inertial detectors. This 
stability implies i) the analyticity and boundedness of 
the Minkowski modes (3) on an inertial trajectory in the 
domain Im(r) < 0. This in turns implies ii) the analytic- 
ity of the Wightman function in Im(r) < 0, iii) A 1 ^ = 
and iv) R 1 ^ = 0. 



These four properties are in fact equivalent. Three of the equivalences can readily be shown directly: i) ii) because 
the measure d 3 k/ujk is Lorentz invariant; i) <4> iii) and ii) <^> iv) by Fourier transform (45) and (54) respectively. One 
can also show directly that the transition amplitudes give back the transition rates, again thanks to the Lorentz 
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invariant measure. 



i) <pg (x a In {r)) analytic & bounded " °4£> /2 " fc ii) Analyticity of Wj„(r) 
FT $ $ FT 

/ = LI of£k/2u, k f i?^< L/ = 

M^Vo lvj \i?-"' L Vo 



The second role of Lorentz symmetry pertains to the 
UA trajectories (6), which are orbits of the generator of 
boosts along the direction z, and are periodic in Im(r). 
The first property implies stationarity, and the second 
property combined with the stability of the vacuum gives 
the thermal spectrum, a) These properties correspond 
respectively to eqs. (17) and (18) of the definition of 
the Unruh modes, the requirement of analyticity fixing 
the ratio of the Bogoliubov coefficients, see eqs (21) and 
(22). Stationarity and stability are also transparent on 
b) the expression (33) of the Wightman function W ua - 
These properties are directly inherited from the Unruh 
(or Rindlcr) modes by Lorentz invariance of the measure 
d 3 /c/|k|, this time expressed in the form of eq. (39). By 
Fourier transform, the periodicity of poles of the Wight- 



man function implies that c) both R± are proportional 
to the Planckian spectrum (exp(27rF/a) — 1) , and the 
stability of the Minkowski vacuum fixes the ratio of the 
rates (the pole +ie contributes to i?_), see eqs. (56) and 
(57). Similarly, as seen on (46), the periodicity of the UA 
trajectory implies that d) oc \A-\, and the stability 
of the vacuum fixes the ratio (this will be clearer in sec. 
II D by comparison with the amplitudes for subluminal 
dispersion relations). Finally, the fact that UA trajecto- 
ries are orbits of the boost generator is responsible for 
e) the proportionality of the Bogoliubov coefficients with 
the transition amplitudes cq. (50), which gives the for- 
mer a physical interpretation. 



Properties a)-d) are also equivalent by Lorentz invariance: from eq. (22), (A7) and (50) we have a) <^ d); a) <^=> 
b) because of the Lorentz invariant measure d 3 fe/|k|; b) c) by Fourier transform; and by Lorentz invariance of the 
measure, one shows [6] that d) =>■ c). 



a) analytic & bounded 
eq.(50) t 



LI of d 3 fc/2w fc 



d) 



aUA,LI i ,UA,LI 



_ e -2-nE/a LI °/ d k/2u> k 



b) Poles ofW u „(T) 
t FT 

c) R U + A - LI /R U _ A ' LI = e -2^/a 



In a nutshell, both roles of Lorentz symmetry are suffi- 
cient to the existence of the Unruh effect. The effect can 
be characterised by either of the four properties a)-d) 
presented in sections I A-I D because they are equivalent 
by Lorentz invariance. We will now show that Lorentz 
symmetry is also necessary for the existence of the Unruh 
effect. 



II. NO UNRUH EFFECT WITHOUT LORENTZ 
INVARIANCE 



Whenever possible, we will establish general results 
valid for arbitrary dispersion relations (DR) which we 



note uJk- The phase and group velocities will be noted 



k 



Vn = 



dk 



(58) 



We recall that we work in the units where the velocity 
of light is 1, which is the asymptotic velocity of the UA 
trajectory (6). The DR is called subluminal (possibly 
on a finite interval only) if v v < 1, and super luminal if 
v v > 1. We will illustrate our results with the particular 
case of linear dispersion relations 



vk , v ^ 1 . 



(59) 



This case can be treated to a large extend analytically be- 
cause the field still enjoys a Lorentz symmetry. The state 
of the field is the ground state defined in the preferred 
frame and noted |0) (the pseudo-Minkowski vacuum). 
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A. Bogoliubov coefficients 

We first extend the construction of I A to linear DR 
(59) and show explicitly that the density matrix p defined 
at equation (16) is not thermal. The sub- and super lumi- 
nal cases must be treated separately. The construction 
of the Unruh modes is identical to the relativistic case 
for supcrluminal DR but it must be amended for v < 1 
because of the instability of the vacuum w.r.t. frames 
(3 > v. We then briefly discuss generalizations to arbi- 
trary DR. 

1. Linear dispersion relations 

The starting point of the analysis of sec. I A is the 
observation that the i?-wedge is globally hyperbolic, or 
equivalently a smooth time coordinate (namely r\ from eq. 
(8)) can be chosen in R such that the surfaces r\ = etc 
are Cauchy surfaces. This is a necessary condition so that 
the solutions of the wave equation Dip = are uniquely 
determined by boundary conditions on that Cauchy sur- 
face. This can be easely adapted to the linear DR (59). 
The wave equation in the preferred frame is 

(d 2 -v 2 (d 2 x + d 2 y +d 2 z ))ip = (60) 

which is trivially put into the form (1) by a rescalling t' = 
vt since the mass-shell relation (59) enjoys an 5*0(1, 3) 
symmetry. The causal properties (60) are determined 
by the pseudo-light cones v\t\ — r, so we define pseudo 
Rindlcr wedges R(L) = {(i, x,y, z)/z ^ 0, v\t\ < ±z}. 
The metric in these globally hyberbolic space-times is 
brought into static form by the introduction of the 
pseudo-Rindler coordinates 

t = - sinh(r)) , z — cosh(?)) (61) 

v 

in R and with a minus sign in L. We can proceed with 
the quantization following sec. I A step by step up to eq. 
(14) and the definition of the a pseudo-Rindler vacuum 

|6*>|6i>. 

At this point we must distinguish between super and 
subluminal dispersion relations. If v > 1, the energy of 
the pseudo-Minkowski modes is in the preferred frame 

^ k = v 2 {k 2 z +^l)>v 2 k 2 z >k 2 z . (62) 

The pseudo-Minkowski modes are therefore analytic and 
bounded functions over the same domain T as the 
Minkowski modes. The pseudo-Minkowski vacuum is sta- 
ble, i.e. this is the lowest energy state in any inertial 
frame. We can therefore continue to follow the procedure 
of sec. I A step-by-step, define Unruh modes, which gives 
the unitary map between the pseudo-Minkowski vacuum 
and pseudo-Rindler Fock state (27). The point of depar- 
ture with sec. I A is that we assume that the observers 
live in the wedge R — {z > 0, \t\ < z}. If the DR is su- 
perluminal, R C R and there exists observables such that 



supp(O^) C R\R where the pseudo-Rindler vacuum |0) ^ 
is not defined, see fig. (1). The trace over L of |0m)(0m| 
is therefore not given by (29). 

If the DR is subluminal, modes with transverse 
wavenumbers small enough to verify 

kl < l —^kl (63) 

break the positivity condition (62). For those modes, 
we must require that ip^ of eqs. (17), restricted on the 
future pseudo-horizon H + , be analytic and bounded in 
Im(y) > 0, where V = vt + z. The unique analytic 
continuation is (V + «e) i0 , which leads to an inversion 
of the ratio $n/&n = e +nQ . (This means not only that 
for observables restricted to this R wedge, the density 
matrix is not thermal, but that it is not a trace op- 
erator since Tr(p) = oo.) Then, as in the super lumi- 
nal case we assume that the observers live in the wedge 
R = {z > 0, \t\ < z) C R and similarly L C L. The state 
resulting from tracing over L cannot be thermal because 
it is still correlated in the region L\L. 

In conclusion, the pseudo-Minkowski vacuum re- 
stricted to the Rindler wedge R is not a thermal state. 
The Bogoliubov coefficients a and /3 give the expression 
of the pseudo-Minkowski vacuum in terms of the pseudo- 
Rindler quanta, whereas we are interested in the expres- 
sion of the pseudo-Minkowski vacuum in terms of the 
Rindler quanta, i.e. 

(0 M \a^a^\0 M ) ^<5(A-A')|/3a| 2 , 
(0 M \a$aj;,\b M )?6(\-\')axP\. (64) 

We will see in sec. II D that the second interpretation of 
the Bogoliubov coefficients in terms of S-matrix elements 
is also lost. 



2. Generalizations 

In the more general case uj 2 = F(k 2 ), the 
mode equation in the preferred frame is given by 
[d 2 - F(d 2 )] ip{x) = 0. Since it lacks an 5*0(1,3) sym- 
metry, the wave equation in Rindler coordinates mixes 
non-linearly time and space derivatives, and the notion 
of Rindler wedges R, L looses its significance. 

But since the settings for the Unruh effect pick a pre- 
ferred direction, see eqn. (6), one could be curious to 
consider instead anisotropic dispersion relations with an 
5*0(1,1) symmetry 

c 2 (fc z ,k ± ) =k 2 z +F(k\). (65) 

This is of course very contrived physically, but we will 
see that formally this mimics rather closely the settings 
of Hawking radiation. The pseudo-Minkowski modes so- 
lutions of (65) are analytic and bounded in T if and only 
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FIG. 1: The Rindler wedges R and L for subluminal (sub- 
script <) and superluminal (subscript >) linear dispersion 
relations (59). 



if F > 0. We can go on defining Rindler coordinates in 
which the wave equation takes the form 



(3* 



dl + e^F(dl)) V 



0. 



(66) 



The solutions of positive Rindler frequency are given by 
(11) with &j_e^ replaced by \J F(kj_)e^ in the argument 
of the Bessel function. We see again that provided F > 
0, the properties of the mode functions are not altered 
w.r.t. the Lorentz invariant case and we can proceed in 
a similar fashion as in sec. (I A). In that case we thus 
recover the Unruh effect at the level of the Bogoliubov 
coefficients, i.e. eqs. (21)-(23) hold. This is because the 
characteristics t = ±z of eq. (66) generate the Rindler 
horizons H ± . We shall return to this in sec. HIE. 



B. Wight man function 

We examine how modifying the dispersion relation af- 
fects the analytical properties of the Wightman function. 
In brief: 1) the analytic properties are essentially the ex- 
pression of the stability or instability of the vacuum in 
all inertial frames. There is therefore a sharp distinction 
between super- and subluminal DR, the latter defining 
pathological models. 2) Both properties of the Wight- 
man function evaluated on UA trajectory (33), that is 
invariance by boost and the equilibrium condition (34) 
are lost. 



1. Inertial frames 

In the preferred frame, the solutions of the wave equa- 
tion ipk oc e - ! "fc t + lkx have a positive Klein-Gordon norm. 



In a boosted frame t' = -f(t+(3z), z' = j(z+(3t), the norm 
of the modes ip k oc exp (— iuj' k t' + ik! z') now depends on 
the sign of the boosted frequency oj' k — 7(0;^ — (3k z ) 



(tp v , ^ fe )cxsgn(u4)<5 (3) (k-k'). 



(67) 



We must thus distinguish two cases. First, with superlu- 
minal dispersion relations 



> 1. 



(68) 



solutions of positive frequency (or norm) in the preferred 
frame have positive frequency in all frames. The vacuum 
|0) defined in the preferred frame is therefore the ground 
state in all the frames. Hence the Fourier transform of the 
Wightman function W = (0\<p(x + y)ip(y)\0) = J flfk 
contains only positive frequencies in any frame, and pro- 
vided the phase velocity does not vanish (otherwise the 
factor \/2uj k introduces poles or branch cuts), W{x) has 
the same analytical properties as in the relativistic case 
in the sense that the regulator e~ ek is equivalent to the 
replacement t n- t — ie. If the DR violates the positiv- 
ity condition (68), the Fourier transform of W contains 
negative frequencies w.r.t. the modes defined in boosted 
frames such that 1 > /3 > v v and the analytic properties 
of W are changed. 

Let us illustrate this with the linear dispersion relations 
(59). The unregularized Wightman function is 



Wz„(r) 



8tt 2 7/^|t| J 



dk 



where we note i] T = sgn(r). The absolute values come 
from r — \z\, see eq. (31). As we said, if v > |/3| we 
can regularize the integrand with e~ ek , or equivalently 
replace r by r — ie. The integration is straightforward, 



1 



1 l-/? 2 

4tt 2 v v 2 — f3 2 r(r — ie) 



(70) 



and one can write the T-dependent part in three equiva- 
lent ways 



1 

t(t — ie) 



1 



1 



(71) 

t z — ir\ T e [t — ley 



The third expression shows explicitely that Wi n is ana- 
lytic in the lower half complex r-plane. For subluminal 
DR and frames such that /? < v < 1, this remains true. 
If v < (3 < 1 on the other hand, the Wightman function 
admits a different representation 



Wm(r) 



8n 2 vjf3\T\ J 



dk ( e - i t k ( v ~l 3r lT)(r+iri T e) 



the integration of which gives 

1 1-/3 2 



Wl n (T) 



4ir 2 v v 2 - j3 2 (t + ie) 2 



(73) 
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It is the complex conjugate of the (70) and is analytic in 
the upper half complex r-plane. As a result, R_ will be 



found to vanish instead of R + . 



2. UA frames 

The same considerations apply to the Wightman function on a UA trajectory. Substituting the parametrization 
(6) into the regularized integral expression of the Wightman function gives 



W ua (S,f) 



87T 2 |Az| 



I fg— i2 j sinh(a<5/2)(u v cosh(ar) — r;a| sinh(ar)|) g — i2 J sinh(ac5/2) (v /p cosh(a-7-)+?j ( 5 1 sinh(ar) | ) \ g — ek ^74^ 

v v I / 



We introduced the notations 

S = n - t 2 , 



n + t 2 



(75) 



and Az = z(ti) — z(t 2 ) = 2 sinh(a(5/2) sinh(af). For 
superluminal DR, that is satisfying (68), the prefactor 
u v cosh(af) ± ris\ sinh(af )| in the first phase is stricktly 
positive, so that e~ cfe can be replaced by S — > S — ie since 

sinh(x — ie) = sinh(a;) cos(e) — i cosh(x) sin(e) 
= sinh(.x) — ie cosh(a;) + 0(e 2 ) 

sinh(a;) - ie (76) 

If the DR is subluminal at some value k, the function 
v v (k) cosh(af) ± 77,5! sinh(af)| changes sign at r(fc) given 

by 



tanh (a\r(k)\) = v v (k) 



(77) 



We cannot replace e~ ek by a single prescription S ± ie. 

The linear DR (59) provides again a good illustration 
of this. Once the integration in (74) done and the frac- 
tions combined, we have 



W ua (S,f) = - 



1 



16tt 2 sinh 2 ( a( 5/2) 
1 



[v 2 cosh 2 (or) — sinh 2 (aT)] — ie 2a At 



(78) 



with aAt = 2 sinh(ac5/2) cosh(ar). If v > 1, we can re- 
place the second denominator by 



[v 2 cosh 2 (aT) — sinh 2 (ot)] x (1 — ie) 



(79) 



and finally absorb the ie into sinh 2 (ai5/2) as done in (76). 
Hence 



W>1, W ua (T 1 ,T 2 )=Wu(S) xf(f) 



(80) 



where the Lorentz invariant Wightman function Wli(^) 
is given at eq. (33) and the function / is given by 



1 



v v 2 cosh (ar) — sinh (ar) 



(81) 



Because of /, W ua is not stationary, and a fortiori not 
thermal. We stress that the factorization of the depen- 
dences in S and f is not generic. It is a consequence of the 
linearity of the dispersion relation, and of the properties 
of the UA trajectories. 

For subluminal DR, the function v 2 cosh 2 (af) — 
sinh 2 (af) becomes negative at the time p given by 
tanh(ap) = v. Thus for f < p we can still use the pre- 
scription 5 — ie, but for f > p we must replace it by 
5 + ie. 



C. Transition rates 

In the Golden rule limit (53), the transition rates are 
given by a Fourier transform of the Wightman function. 
The properties of the latter therefore pass on to the for- 
mer: the transition rates are not stationary, and their 
ratio is not the Boltzmann factor. We show this explic- 
itly with linear DR (59). In appendix B we use asymp- 
totic expansions of the Wightman function and transition 
rates to estimate the proper time interval elapsed since 
the beginning of the acceleration after which stationarity 
and thermality are lost for more general DR. The result 
is 



2tt < at < - In 



M M\ 



(82) 



The lower bound excludes transients. If we take M equal 
to the Planck mass, the upper bound can be as high as 
100. 



1. Inertial detector 

We start with the integral representation (69) regular- 
ized with e~ ek and integrate by parts 
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• 1 pOO 

Wm(r) = [G(fc) (e— - e-^) e'^ + ^ J Q * G ( fc ) ~ /^"^ - K + /3)e-^ +T } e" efc (83) 

where G(k) = f k f^ and we note cu± = j(cuk ± /3k). The boundary term vanishes at both the lower and upper bounds. 
We then substitute this expression in (53) and exchange the order of integrations 

/+oo 
dr e* iET W In (r) 
-oo 

= \ dkG{k){(v g -p)5{u-±E)-(v g + /3)6(u + ±E)} e~ ek , (84) 



As expected, R+ 1 ^ if cj_ < 0. We can further write 

Rl - = £z E ( s p G (V) - 4C(k+)) (85) 
p 

where fc^ and are defined by 

£ = w± (fc± ) , s± = sgn(« 3 (fc± ) ± /?) (86) 

The calculation with a linear dispersion relation (59) 
provides an independent check of this result since in that 
case we have integrated over k first, and the integral over 
t can be done by application of the theorem of residues. 
For superluminal DR, we use (70) and get 

R+ — , R- = ^- ^ , (87) 

2-7T V{V Z - (} Z ) 

which is also the result obtained by application of the 
expression (85). 

For subluminal DR, the previous result is still true as 
long a j3 < v. For /? = v the Wightman function is not 
defined and for (3 > v the ie prescription must be complex 
conjugated, i.e. the numerator of (32) is (v 2 — /3 2 )(r + «e) 2 
and we now find 

R -=°> R+ = 9 ^ v (^ v 2-y (88) 

which matches with (85). The vacuum of the quantiza- 
tion in the preferred frame appears to this observer as a 
negative energy state with respect to his ground state. 



I 

2. UA trajectories 

We lack analytical tools to study the general case of 
an arbitrary dispersion relation. In appendix B we use 
Taylor expansions to show that stationarity is lost after 
a few thermal periods a/2ir at best. Below we estab- 
lish results for the linear dispersion relations (59). The 
Wightman function (74) is not stationary, which means 
that we cannot calculate the transition rates (53) with a 
contour integral. We find however that the "even" part 
of the transition rates 

R±(r) = 1{R±(t)+R±(-t)} 

/CO 
drie^^T-n) (89) 
-oo 

can still be calculated as a sum of residues. We call it 
the "mean rate". If R is not constant, then dR/dr is not 
an even function and R(t) is not constant. 

We must as before distinguish between super and sub- 
liminal DR. We consider first phase velocities v > 1. We 
introduce the Rindler time ap defined by 

tanh(ap) = v (90) 

in order to write the denominator in (81) as a prod- 
uct u 2 cosh 2 (af) — sinh 2 (af) = (v 2 — 1) cosha(p — 
f) cosha(p + f). After the change of variable x = 
the expression of R + is 



R = g 2 asinh 2 M Rc [+°° ^ 



/+oo P t T x 1 
dx (91) 
-oo sinh 2 (x - ie) cosh [x - a(r - p)} cosh [x - a(r + p)] 



+ 47r 2 2v J_ 00 sinh 2 (x - ie) cosh [x - a(r - p)] cosh [x - a(r + p)) 

Wli has a familly of douple poles at inn + ie and / has two families of simple poles 



x p = a(r - p) + + ipw , y p = a(r + p) + i^ + ipir (92) 



12 



Note that these poles are shifted by with respect to the double poles inn. The importance of this structure will 
be better understood in the calculation of the amplitudes. One obtains 

R±(T) = RL±If[T) + /* (M2E(r + p)] M2E(r- P )}\ 

W ± y> 16tt 2 v 2 sinh (irE/a) \ cosh 2 a(r + p) cosh 2 a(r + p) J V ; 



where /(r) is defined at eq. (81). We refer to the figures 
2 and 3 for further details on these expressions. 

Although R±(t) is not the transition rate, the fact 
that the second term (93) is common to both mean rates 
shows that the rates are not time independent and that 
their ratio is not a Boltzmann factor for times |r| > p 



R+ , ( 2itE 

— i= cxp 

R- \ a 



(94) 



We also calculated numerically the transition rates. They 
are shown on fig. 4 and 5. We observe that they are 
equal to the Lorentz invariant ones inside the interval 
|t| < p, and asymptote rapidly to zero outside, but are 
equal. Between these two regimes, the transition rates 
experience a burst. 
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FIG. 3: The ratio (R+(r)/R- (r))e 27ri5/[! of the mean rates 
normalized to the Boltmann factor. We took v = 1 + 
10~ 10 ,a = 1, and E = 1 (plain), E = 0.8 (dashed). Devi- 
ations from the Lorentz invariant result are of order 1. 




FIG. 2: Dependence of R+(t)/R+ on the velocity. The two 
curves correspond to v = 1 + 10 -10 (plain) and v = 1 + 1CP 8 
(dashed) with a — 1, and E = 1. The deviations from the 
Lorentz invariant occurs at times |r| > p where p is defined 
at eq. (90). Its values are respectively ~ 11.9 and ~ 9.6 



For v < 1, the Wightman function does not have a sin- 
gle analytic expression in S for all values of t and we can- 
not calculate R± by an integral contour. In preparation 
of the analysis of the S'-matrix elements, it is nevertheless 
useful to examine the changes in the analytic structure 
of the integrand (89) compared to the superluminal case. 
For v < 1, we introduce in place of (90) 



coth(ap) = v . 



(95) 



The rapidity ap was previously introduced after eq. (81) 
as the critical value of proper time after which the pre- 
scription in the Wightman function must be changed to 



FIG. 4: The transition rate R+(t) normalized to the Lorentz 
invariant value. We took v = 1 + 10 _6 ,a = l,a = 0, and 
E = 0.1. It is equal to 1 for times |r| < p — 7.3, and after a 
burst around \r\ — p it vanishes exponentially. 



r) sinh a (p + r), whose poles are 

x = ci(t — p) + ipn 



Vv 



a(r + p) + ipir 



(96) 



ic 



The denominator of / is now (1 



i sinh a(p • 



Contrary to superluminal velocities, all the poles are ar- 
ranged on lines Im(r) = ipn. We will now see how this 
analytic structure appears in the transition amplitudes. 



D. S'-matrix elements 

We first give exact expressions. Again sub- and super- 
luminal DR give different results. The ratio of the transi- 
tion probabilities is the Boltzmann factor if Uk > k, and 
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FIG. 5: The ratio R+(t)/R-(t). The dashed line is the Bolt- 
mann factor. We took v = 1 + 10~ 6 ,a = 1, and E — 0.1. 
Interestingly the rates are equal for times |r| > p. 



one if ujk < k. The latter result is a dramatic consequence 
of the instability of the vacuum. 

We then calculate again the transition amplitudes by 
the method of steepest descent for the role of the analytic 
properties of the modes in the previous results appears 
more clearly with this method. 



1. Exact expressions 

The inertial case is readily dealt with. The amplitudes 
are still given by (42). Provided condition (43) is verified, 
A 1 ™ vanishes and A 1 ! 1 ^ for the wavenumbers allowed 
by the conservation of energy. If v v < 1 over some inter- 
val, the energy of the scalar quantum in the rest frame 
of the detector can become negative, thus opening the 
channel A+ l ^ while closing the other A 1 ! 1 = 0, which 
is again interpreted as the instability of the vacuum for 
frames such that (3 > v v . 

We now turn to the UA trajectory. We use the same 
notations as in sec. IC except for k± in eq. (46) 
which should be replaced by its more general expression 
z = ^/k + k^, see eq. (44). The super luminal case is a 
repetition of the Lorentz invariant one because uju > k z 
in every Lorentz frame. We get 



a ua = ~ i2 9 e T 21 
±M ^2u k {2nf a 



In -\f— ~k~Jk~^ we now have 



A 



UA 
±,k 



oc 



±iE/2 



dye 



±iEy—iz cosh(y) 



(99) 



where we note z = ^—k + k-. This integral is 
defined as the analytic continuation of the modified 
Bessel function in terms of Hankel functions, K v (z) — 



-i\e- l ™l 2 U v (e"/ 2 z) 



aua _ 

•^i.k — 



-gn 



g-t- 2 



V^k(2ir) 3 a 



H™ M (z). (100) 



Since H_l(z) = e~™"Hl (z), we now have 



pUA 
pUA 



= 1 , 



(101) 



This result could have been almost anticipated from the 
response of the inertial detector since then we also no- 
ticed that both channels, excitation and desexcitation of 
the detector by emission of a quantum, can occur. What 
is perhaps surprising about (101) is that it is indepen- 
dent of the energy gap, as if the detector was coupled to 
a "reservoir" of infinite energy. We explain this by the 
fact that the vacuum defined in the preferred frame is 
unstable in frames /? > tanh(w ¥ ,), hence along most of 
the UA trajectory. (The ratio (101) is what one would 
obtain in a thermal bath at infinite temperature, but one 
should not use this misleading analogy because the UA 
detectors does not react as in a thermal bath.) 



2. Steepest descent approximation 

Looking at the integral expressions (46) and (99), we 
see that they differ only by a shift in/ 2 of the variable 
of integration, since isinh(y) = cosh(y + in/2). A simi- 
lar shift by i7r/2 was found in the position of the simple 
poles of the Wightman function between the superlumi- 
nal (92) and subluminal DR (96). To better understand 
how these two analytic structures are related, we evaluate 
the previous integrals by the method of steepest descent. 



K 



■XsfKkZ). (97) 



dx 



AE/a)f(x) 



/2tt 



The only essential difference is that fc + fc_ = u>% — k 2 z ^ 
k\. But since K_ v (z) — K„(z), the ratio of the proba- 
bilities is still a Boltzman factor 



fix) 



irjx 



(JL 



(102) 



pUA 

+ ,k _ -2TrE/a 
pUA e 

-,k 



(98) 



The description of the subluminal case with small 
values of the transverse wavenumbcr is however differ- 
ent. Now either k + or fc_ is negative. Let us say 
/c_ < for defmiteness. With the change of variables 



where r/ = ±1. From now on we work in the units of 
E. The results of the analysis, presented in appendix 
C, are summarized by the fig. 6 and 7. The two cases 
they represent are respectively A > and < 0, where 
A = 1 + k 2 z — u) 2 . The figures represent the real part of 
f(x + iy), the dots are the saddle points x± (two in each 
case), and the curves are the path of stationary phase, i.e. 
the solutions of Im(/[7(A)]) = lm(f(x±)). Only in the 
case A < does this curve verify the additional condition 
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arg (f"(x±)(x\y — x±) 2 ) = 7r so that it corresponds to 
the steepest descent path. 

In more details, for A > (subluminal DR and slightly 
super luminal DR), the steepest descent path does not 
exist. The curve of constant phase first climbs from an 
abyss before reaching the saddle point, and leaves it to 
ascend a mountain. If it existed, the steepest descent 
path would instead continue progressing toward increas- 
ing Im(x) in order to reach the abyss nearby. 




FIG. 6: Region I and II of parameter space (lo, k z ) corre- 
sponding to A < 0. The real part of f(x) as a function of 
Re(;r) and Im(x). The dots are the saddle points and the 
curves on the surface are the solutions of (C7). We took 
uj — 1.5, k z = 2 and ij — — 1. The figure for ij — +1 is similar. 

For A < (superluminal DR), the steepest descent 
path exists. It is the curve of stationary phase (C14) 
passing through the saddle point x + . The crutial differ- 
ence with the previous case is that abysses on the side 
He(x) > face abysses on the opposite side, while in pa- 
rameter regions I and II their relative positions is shifted 
by 7r/2, mountains thus facing abysses. The real part of 
f(x + iy) is indeed 

fc_ < , Re(/) = —2772/ + C smn f x + - m -/—) sm 2/ > 

V 2 k + J 

k- > , Re(/ ) = — 2r/y — C cosh ^x + ^ In ^— ^ sin y . 

(103) 

The behaviour at large x of the former expression, valid 
in regions I and II, depends on the sign of x. For x > 
0, it diverges to +oo for y e ]ir, 2w[, and for x < it 
diverges to +oo for y e ]0, ir[. Again, mountains arc 
facing abysses. The asymptotic behaviour of the second 
line of (103), relevant for region III, is independent of the 
sign of x and mountains on one side thus face mountains 
on the other side. We note finally that both the relative 
positions of the mountains and abysses and the angle of 




FIG. 7: Region III of parameter space corresponding to A > 
0. The real part of the phase f(x) as a function of Re(;r) 
and Im(a;). The dots are the saddle points and the curves on 
the surface are the steepest descent paths (plain for x + and 
dashed for We took w = 7, k z — 5 and r\ = —1. The 

figure for r) — +1 is similar. 



the tangent at the saddle point go hand in hand since 
f(x) = irjx + f"(x), and it is /" which determines both. 

This establishes that the additional poles of the Wight- 
man function (92) and (96) are in one-to-one correspon- 
dance with the positions of the maxima of Re(/) (again, 
this was expected since transition amplitudes and Wight- 
man function are build from the modes). 

We finish with the expression of the amplitude in re- 
gion III evaluated on the steepest descent path 



-(£/a)(7r-arctan y/-A)-(E/a)\ 



iq=-l — 



fcf)V4 



iE/2a 



(104) 



e 2 for z > 1, these exprcs- 



T i - e * E / a T 

Reminding that K v ^ 
sions are indeed the uo 3> k z limit of (97). 



E. Discussion 

We saw in sec. I that Lorentz symmetry assumes a 
double role: it ensures the stability of the vacuum in any 
frame and endowes the UA trajectories with special prop- 
erties which gave in turn the stationarity and thermality 
of the Wightman function and transition amplitudes and 
rates. In this section we first saw that stability is pre- 
served only by superluminal DR (68). Second, station- 
arity and thermality are lost simultaneously. Third, the 
transitions are equal to the Lorentz invariant ones over 
a finite interval and deviate significantly from them out- 
side that interval. This interval is fixed by the value of 
the Lorentz factor (7) such that the detector probes the 
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non trivial properties of the dispersion relation. Quite re- 
markably, no matter how small the deviation from u = k 
may be, once T is large enough the transition rates de- 
viate appreciably from the Lorentz invariant ones. With 
a linear DR (59), fig. 4 and 5 show that the transi- 
tion rates are exponentially close to the Lorentz invari- 
ant ones for |r| < p and exponentially small and equal 
for |t| > p. The transition between the two regimes is 
moreover sharp. For more general DR, we showed in 
the appendix B that the corrections are controled by the 
Lorentz factor times E/M or a/M, see eqn. (B9). 

By Lorentz symmetry we have four equivalent ways to 
demonstrate the Unruh effect. This equivalence is lost 
in the present case. In particular the mathematical con- 
struction which consists in the definition of a quantum 
field theory in a Rindlcr wedge has no operational mean- 
ing (the Bogoliubov coefficients are not proportional to 
the transition amplitudes of detectors). Note also that 
one should not interprete the value of the ratio of the 
transition amplitudes for supcrluminal DR (98) as the 
proof of the Unruh effect, because only inclusive proba- 
bilities are measurable. In relativistic theories we can use 
the ratio to characterize the effect because, by Lorentz in- 
variance of the measure d k/u>k, the amplitudes give the 
transition rates (56). Without Lorentz invariance, the 
transition amplitudes sum up instead to non stationary 
transition rates. 

In conclusion, we have seen that Lorentz invariance is 
both necessary and sufficient to the Unruh effect which 
can be characterized in four equivalent ways, although 
only the transition rates of a detector are measurable 
quantities. Without Lorentz symmetry, this equivalence 
is lost and the transition rates differ significantly from 
the Lorentz invariant ones after a proper time at which 
the Lorentz factor of the detector is high enough to probe 
the non trivial features of the dispersion relation (v — 1 
or E/M). One expects however that this value is very 
high, simply because the scale at which Lorentz invari- 
ance might occur is very high. It seems therefore unlikely 
that the Unruh effect could be used as a practicle test of 
Lorentz invariance. 



III. COMPARAISON OF THE UNRUH EFFECT 
WITH HAWKING RADIATION 

Hawking radiation (HR) designates the property of the 
vacuum as defined by observers in free fall near the hori- 
zon to correspond to a thermal bath for distant observers. 
The first four subsections contain review material. Their 
presentation is however sketchy and the reader will find 
further details in references [2, 11, 13, 14] in particular. 
The comparison with the Unruh effect is done in section 
HIE. 

To begin, it is worth recalling that HR has little to 
do with gravity in the sense that whether the metric is 
a solution of the Einstein's equations or not is irrelevant 
to the matter [11] (but field equations are necessary to 



establish the laws of black hole thermodynamics). What 
does this mean exactly for a Schwarzschild black hole? 
The metric outside the black hole in static form is 



ds 2 



-f(r)dt 2 



dr 2 
W) 



f(r) = 1-f • 



(105) 



The curvature is given by R ab = 0, C abcd C a b c d = 
12r 2 /r 6 , and the surface gravity is k = GM/r 2 . Hawk- 
ing's result uses only those properties of space-time which 
survive in the limit 



2GM 



cte , M -> oo , G -> 



(106) 



where gravity is decoupled while the geometry is fixed. 
Indeed in this limit the curvature and surface gravity, 
which are defined and calculated without recourse to the 
field equations, depend only on r s and are therefore in- 
variant. Two additional quantities are important (the 
field equations are required to show the following results) . 
First the evaporation time t cv cx G 2 M 3 /he 5 — > oo, that 
is backreaction vanishes, which shows that the limit of 
fixed background is self-consistent. Second, the entropy 
S = AttGM 2 /he — > oo, hence quantum field theory is in 
principle valid down to arbitrarily small scales. This is 
because to a region of size r s is associated a maximal 
density of state e s [12] in a dynamical theory of gravity. 
This statement is incompatible with a relativistic field 
theory which has an infinite density of states by Lorentz 
invariance. Conversely if relativistic quantum field the- 
ories in curved space time are valid, the entropy of the 
black hole (which saturates Bekenstein's bound) must be 
infinite. 



A. Radially free falling frames 

Let us now proceed with the derivation of HR. As 
we said the state of the field is the vacuum as seen by 
geodesic observers. Begin thus with the construction of a 
one-parameter family of coordinates (r, p, 9, <f)) attached 
to observers radially free-falling [13]. This family is pa- 
rameterized by the Lorentz factor 7 of the observers at 
infinity 



7 2 = 



1 



1 

P 



(107) 



r is their proper time and p is the proper radial coordi- 
nate on the surfaces of constant r. The radial velocity is 
given by 



v(p) = -y/l-pf[r(p)]. 



(108) 



The sign in front of the square root is negative be- 
cause we consider infalling geodesies. These coordi- 
nates are related to the Schwarzschild coordinates by 
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dr = 7 (dt — jdrj and dp = dr/j, or in integral form 

T = 1 {t + R{r)) (109) 



with 1 
R(r) 



P = 

1 

rv{r) + r s In 
(P-2) 



r„ In 



P+2- (l-p-y/l-pvj 

To \ / 



r , In 



p/(r) + 21- 



+ etc 



(110) 



Only the near and far horizon limits are actually inter- 
esting 



R{r)=r s \n[ r -^\+C + 0" 



r^oo, J2(r) = + C + O (y) . (Ill) 

Near the horizon the R(r) ~ r*(r) = r + r s ln(r/r s — 1), 
the tortoise coordinate defined by dr* = dr/f(r), so 
r — > 7(t + r*) where i + is the advanced Eddington- 
Finkclstcin null coordinate. Far from the horizon, r — > 
7(* — Woo?') is naturally the proper time of observers 
boosted with the static observers w.r.t. the black hole. 
The line element 

ds 2 = - (l-v 2 ) dr 2 -2vdrdp + dp 2 +r 2 dfl 2 (112) 

is stationary and rotationaly invariant. The radial com- 
ponent of the shift vector is —v(p). The Schwarzschild 
radius r — r s corresponds to v(p s ) = —1 and near the 
horizon the expansion 

v = — 1 + x + 0(x 2 ) , x = K 7 (p — p s ) (113) 

will be used, where 

*r = ^ (H4) 
is the surface gravity seen by these observers at infinity. 

B. Partial wave decomposition of the field equation 

We now pick one of these free fall frames and define the 
modified dispersion relation from the quadratic action 2 



1 



dr d 3 x^{(n a d a ^ 



written in covariant form with the help of the 1 + 3- 
decomposition of the metric: n = (l,v,0, 0) is the 
unit vector normal to the surfaces of constant r, q ab = 
g ab + n a n b is the induced (contravariant) metric tensor 
on these surfaces, and A = q~ x l 2 d a (yfqq ab db) the corre- 
sponding Laplacian. In the following two subsections we 
review in turn the Lorentz invariant case F = 0, and the 
changes introduced by F ^ 0. 
The Klein-Gordon scalar product 

(&V0=* Jd s x^q{4>*{n a d a ^)-^{n a d a <l>)*} (116) 

is conserved on the solutions of the field equation. This 
equation 

d b (^n b n a d a ^) = A0 + F(A)(p (117) 
is separable for the ansatz 



olra 



-Yi m (6,4>) <p u i{p) 



(118) 



where the radial functions tp u t are solutions of 

(v 2 - l)ip" + 2v(v' - iu)tp' - lj(uj + iv')<p 
= [Vt+FiD)}^. 



(119) 



The prime stands for the derivative with respect to p. 
The auxiliary functions appearing in this equation are 
the transverse momentum 



A: 



2 1(1 + 1) 



(120) 



and the effective potential V; = kf — 2vv /r ^. We intro- 
duced the second order derivative operator D defined by 

A^ (m = r~ x Y lm {D<p ul ) (121) 

Its explicite expression is Dip^i = (d 2 — kf) <p u i. 



HR with Lorentz invariance 



In this section, F = 0. Equation (119) has the partic- 
ular property that the coefficient of d 2 ip vanishes at the 
horizon. It must therefore be solved separately on each 
side of the horizon. We thus expect two classes of solu- 
tions, the ones that are singular at the horizon and the 
ones which arc regular. They will describe respectively 



+ q ab d a ct>d b 4> + 0F(A)0} (115) 



1 We report a typo in eq. (3.4) of [13]. One should read r/7 = 
T = t + rVl -pf + ■■• instead of T = t + r(l - pf) + .... It 
can be easely checked in taking the limit p — > 00 that the latter 
expression is incorrect. 

2 It is customary to choose the frame at rest at infinity as the pre- 
ferred frame. This is probably an implicit simplifying assump- 



tion, or the natural choice in the idealized model of an isolated 
black hole. There is otherwise no reason to opt for this frame, to 
which one could prefer for instance the rest frame of the cosmic 
microwave background. We shall therefore not follow the cus- 
tum. Besides, this does not introduce any complication because 
only the asymptotic behaviours (111) of the free fall coordinates 
and modes (126) and (128) arc relevant, and they are universal 
up to a trivial factor of 7. 
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radiation at infinity and the free fall vacuum. HR is es- 
tablished by showing that these two sets of modes are 
related by a Bogoliubov transformation similar to (21)- 
(23). 

This is most simply done by solving the field equation 
in the eikonal approximation 

0= -A(T,p,d)e tS ^ e \ g ab d a Sd b S = (122) 
r 

where A is a slowly varying function compared to the 
Hamilton- Jacobi action S. By separation of variables, 
S = -ujt + + 1)0 + J p dp'k u (p') where fc u are the 
solutions of the quadratic equation 

(w - v{p)k u f =kl + kf . (123) 

with ki defined at cq. (120). The two roots are 

k±(u,p) 



l-v 2 (p) 



-ujv ± ,Juj 2 + (v 2 - l)kf^ (124) 



Far from the horizon v — > v cai ki — > 0, and k± — > to/ (uco± 
1). The solutions fc_(c<j) < are propagating towards the 
black hole (thus describe infalling radiation) 



p -» oo , mf -» — cxp { -iu) \t H — 

r { V 1 - Woo 

= "^L e -iju(t+r) 

r 



(125) 



To get the second expression we substituted eq. (111). 
The solutions k + (uj) > propagate away from the black 
hole and describe the HR emerging at infinity 



p — > oo . 



/hr 



A 

— exp < — IbJ I T 

r I V 1 + Woo 



A 



-iju}(t—r) 



Near the horizon, the infalling solution becomes k- 



I + & + 0(x) and 



P^ P~t , 



^inf _^ e -iw(r+p/2) 



(127) 



and does not present any particular interest. The state 
of these modes is assumed to be the vacuum. The other 
solution k + <~ ^ is singular at the horizon and describes 
two types of modes with support on either side of the 
horizon 



(128a) 



-hr 



— > 9(x) exp I —iuj It — In a; 



= 9(r -r s )e- i ^ t - r '\ (128b) 
iP t„* ^ cxp ^ r _ J_ ^(-^ | 

= 0(r s -r)e-^( 4 - r ^ (128c) 



where we have omitted terms 0(1) in the phase. The 
modes <fr hr describe Hawking radiation. Note that the 
asymptotic forms (128b) and (126) can be written e~ %1U)u 
in terms of the advanced null coordinate u = t — r*. 
The modes <fi ptn * are trapped inside the horizon and de- 
scribe the partners of the quanta of Hawking radiation. 
They are complex conjugated so that they have a positive 
norm. 

This completes the description of the solutions rele- 
vant for the quantization by observers at infinity. Freely 
falling observers on the other hand are assumed to expe- 
rience the vacuum as they cross the horizon. One defines 
the free fall vacuum in a similar way to the Unruh vac- 
uum in sec. I A. One constructs a basis of modes <j> u 
regular across the horizon and with the following prop- 
erties: 

i) they are eigenmodes of the Killing vector d T , 



J uj\m 



-Y lm {6,4>)<P U M 



(129) 



ii) they are eigenfunction of the Lie derivative with re- 
spect to the unit vector n orthogonal to the surfaces of 
constant r 

£n €lm = n a d a ^ lm = -ifi(p) 4% lm ■ (130) 

with £1 > 0, so that they have positive norm (116). One 
notices that near the horizon we have the identity 



£, n 6 



-i-yiAj(t— r*) 



-iju>(t—r*) 



(131) 



Hence the solutions e _I7 "(* _r *) have positive (negative) 
free fall frequency outside (inside) the horizon. Similarly 
the solutions e + l 7"(*-»"*) have positive free fall frequency 



(126) inside the horizon. Reminding that e 1 ' 



'< one 



infers that near the horizon 



owe-*"* \x\ iu '^ 



0(x) + ^9(-x)\ (132) 



(we suppressed the term Yj TO for a better lisibility). The 
unique analytic continuation of x %u / K ~< which is bounded 
in the domain {Im(i) < 0, Im(r,) < — Im(t)} is 



<P%l = a^i(x + ie) 



(133) 



The coefficient a is fixed by the normalization. With the 
branch cut of the logarithm along the negative real axis, 
the r.h.s. evaluates to the sum 



C/ m = «^0 hr (p) + ^0 ptn *(p) 

with the ratio of the Bogoliubov coefficients 



-27zuj/k 1 



(134) 



(135) 
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D. HR without Lorentz invariance 

There are three key elements in the derivation of (134). 
One is the universal behaviour of the coordinate r(t,r) 
near the horizon. Modifying the dispersion relation does 
not change this. The second is the logarithmic depen- 
dence of the modes near the horizon as a function of x. 
The third is the branch cut associated with the analytic 
extension of this logarithm. It was introduced when we 
chose the state to be the free fall vacuum. Let us now 
examine the changes caused by F ^ 0. The following 
discussion is mainly qualitative. 

The origin of the logarithm is the " kinetic" term of eq. 
(119), i.e. the differential operator on the l.h.s of this 
equation. We can factorize the latter into an infalling 
and an outgoing part di n {dh T <p where 

dinf — d T + d p v — dp — > d T — 2d x (136a) 
<9hr — d T + vd p + dp — > d T + nxd x (136b) 

The solution of the former is <p lni of eq. (127), and the 
solution of the second describes Hawking radiation (128). 
Since this operator is the difference of the l.h.s. of (117) 
with the Laplacian A(f>, whenever one replaces in the ac- 
tion 4>A(j) by cj>(A + F(A))<j> as in (115), this kinetic op- 
erator is preserved. So we also expect to find solutions 
with a branch cut in that case. 

Indeed, taking the Fourier transform w.r.t. x of the 
limiting expressions of (136) near the horizon gives 

i(w + 2k) [nd k + {ito + n)]<p= [Vi + T{k 2 , k 2 )} (p (137) 

where Vi is a constant near the horizon. Note that 
the function T differs from F of the dispersion relation 
by terms containing derivatives of kf, possibly coupled 
to derivatives of (p. For instance, a term D 2 tp gives 
[k 4 + 2kfk 2 + (kf - d 2 kf)] ip. Similarly, a term D n ip 
produces a homogeneous polynome of order 2n multiply- 
ing (p. Only for the s-wave (or in 1 + 1 dimensional mod- 
els) do we have ^(k 2 ,^) = F(k 2 ). These terms could 
affect significantly the grey body factor. Let us now con- 
sider the left hand side of (137). In relativistic theories, 
we know that first black holes are black bodies, so the 
typical frequency of Hawking radiation is given by the 
temperature, i.e. to — O(k), and second, the physically 
interesting region is n^x <C 1. In other words HR corre- 
sponds to low frequencies and high wavenumbers k lo. 
We thus replace u> + 2k on the l.h.s. of (137) by 2k and 
solve the equation with the ansatz (p = (fo(k) x(fc), with 
<p = e^k-^/^- 1 solution of [nd k + (ito + k)] ^ = 
defining the free fall vacuum (since <po(x) °c (x+ie) 1 ^' ' K ' 1 ), 
and x solution of dlnx = —iTj2nk 2 . This approxima- 
tion, also adopted in [14], amounts to neglect the coupling 
between the outgoing and infalling solutions (the latter 
corresponding to the root to + 2k = as we know from 
the WKB solution). This coupling can indeed be argued 
to be innocuous [15]. Under certain assumptions, e.g. 
analyticity of F, the inverse Fourier transform can be es- 
timated in the steepest descent approximation, in which 



case HR is found, see [14] and [16] for more details. Its 
origin is clearly identified as the branch cut of ip . 

This leaves the question of the state. At least one con- 
dition seems necessary so we can assume that the field 
is in the free fall vacuum, namely that the evolution of 
the modes is adiabatic [2, 14]. This places certain con- 
strains on the dispersion relation. For instance, if F is 
polynomial of order 2n, the modified dispersion relation 
(to — vk) 2 = k 2 + kf + F(k 2 ,kf) posesses 2n solutions 
amongst which 2p > 2 are real and 2(n — p) are complex 
conjuguate. In that case a necessary condition for adia- 
baticity is the absence of level crossing between the real 
roots. This can happen via a kind of seesaw mechanism 
if the modifications of the dispersion relation arc charac- 
terized by a very high scale M > w ~ k. It requires 
some care to analyse the contribution of the complex 
roots, but again they should not affect the low energy 
part of the spectrum, adiabaticity implying their decou- 
pling from the high energy modes. A detailed analysis 
of these complex roots for the DR k 2 ± k 2 /M 2 confirms 
this qualitative argument [17]. 

The fundamental part played by adiabaticity should 
not be surprising if one recalls that modified actions such 
as (115) describe an effective field theory, and in a non 
trivial background, both scale separation and adiabatic- 
ity are necessary to validate this framework [18]. 



E. Comparison 

Let us finally return to the question of the relationship 
between the Unruh effect and Hawking radiation. HR 
is habitually proved by finding the ratio of the Bogoli- 
ubov coefficients as in eq. (135) and is therefore deter- 
mined by the solutions of the field equation. One set of 
modes (and therefore the corresponding state) are regu- 
lar across the horizon (Minkowski/Unruh-likc), while the 
other set (Rindler/out going) have a logarithmic singular- 
ity (in V = t + z or x). 

The essential difference is the role played by the hori- 
zon in the dynamics of the field. There are actually two 
notions of "horizon" that should be distinguished. One 
is the surface of infinite redshift associated with the ob- 
servers, that is the surface v = — 1 for observers far from 
a black holes, as examplified by (128) and (131), and 
the null planes t — ±z for the uniformaly accelerated ob- 
servers. The other is the locus of the logarithmic singu- 
larity of the modes. Without Lorentz invariance, these 
two notions still coincide for a black hole, but they differ 
in flat space. Indeed, whether the dispersion relation is 
relativistic or not, the field equation (119), or (136b), is 
singular at v = — 1. In contrast, the horizon of unifor- 
maly accelerated observers is not the locus of singularity 
of the (pseudo-)Rindlcr modes, which is v$t — ±z for 
linear dispersion relations ujk — v^k, and which is not de- 
fined for general dispersion relations because of the non 
linear mixing between Rindlcr coordinates. This explains 
also why the dispersion relations (65) with a SO (1, 1) 
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symmetry mimic better the black hole context, because 
in that case the modes are still singular on the observer's 
horizon. 



As any phenomenon of pair creation from an unsta- 
ble ground state, Hawking radiation is characterized by 
a branch cut, namely (133), or more generally fc _KJ / K -v 
from (137). This branch cut is robust because the higher 
derivative terms do not mix with the operator on the 
left hand side of (119) or its Fourier transform (137). 
We already gave a necessary condition for this: that the 
gradients term in the action be replaced by A + F(A). 
There is a second condition that we did not mention so 
far, although it should be quite obvious. All the results 
of the previous sections depend on the fact that the ob- 
servers at infinity and the observers near the horizon be- 
long to the same referential, that is they are all freely 
falling observers characterized by the same Lorentz fac- 
tor 7. This should be contrasted with the Unruh effect 
where it is necessary to boost a detector continuously (at 
a constant acceleration). In the case of HR on the other 
hand, the redshift between free fall observers near and far 
from the horizon is purely gravitational, that is caused 
by the curvature. An observer equiped with a two-level 
detector and freely falling but with a different Lorentz 
factor 7' would observe similar phenomena as the ones 
described in the first part of the paper. For instance, if 
the dispertion relation is subluminal, it would perceive 
the free fall vacuum of the other observers as unstable if 
it is sufficiently boosted w.r.t. them. 



One sometimes invokes the equivalence principle as 
the reason for the analogy between the Unruh effect and 
Hawking radiation. This is clearly not correct, for other- 
wise we would not expect to find a Hawking-like radia- 
tion of phonons in a variety of condensed matter systems 
(dumb holes) with an acoustic horizon, since this predic- 
tion does not require either the Einstein's equations, nor 
even Lorentz symmetry. What the equivalence principle 
does imply however is that, if the preferred frame is not 
the one with 7 = 1 (and there is no reason why it should 
be, see footnote 2), a static observer at fixed radial dis- 
tance from the black hole should record transition rates 
similar to those of sec. II C. 



In brief, the prediction of Hawking radiation rests on 
the fact that the field equations in a black hole metric are 
singular on a surface which coincides with the horizon of 
asymptotic observers. As long as modified dispersion re- 
lations do not alter this property of the field equation 
(and provided the evolution of the state is adiabatic), 
Hawking radiation is expected to be robust (in the pre- 
ferred frame). By contrast in 3+ 1 Minkowski space, the 
locus where the the Rindler modes are singular coincides 
with the horizon of uniformaly accelerated observers only 
if Lorentz invariance is assumed. 



Appendix A: Bogoliubov coefficients 

Since Rindler and Minkowski modes are solutions of 
a linear equation and form complete families, they are 
related by linear transformations 



dk z (« + /3 A \<*) 

-OO 

/ + OO 
dk z (a L Xkz <p™+Px kz <*) ( A1 ) 
-oo 



From the property (13) one deduces 



4 2 



(A2) 



The Bogoliubov coefficients a and fi are given by the 
Klein-Gordon products a R = (<p M , ip R ) and f3 R — 
((p M * , ip R ) . Let us present the calculation of a R k ^ . Since 
the Klein-Gordon product is independ of the time argu- 
ment of solutions of the wave equation, we evaluate it 
at on the horizon H + . After an integration by parts we 
have 



<*aL 



i2 



dV^dvtf + bnd 



= 2AACAA 







3 dV 



e iu - v / 2 V~ iX +bnd. (A3) 



We note Af the normalization constants of the modes. 
The boundary term is [ ( Pu* ( P\] v Z ■ To obtain the 
second line we used ip M = A/" w e~ iw - y / 2 on % + = 
{U = 0, V > 0} with k- = to k - k z (see eq. (44)), and 
the second term in the asymptotic expansion (20) (we 
assume again narrow wave packets in A). The integral 
can be written as the limiting value of Euler's T function 

f°°dt 

- z T{z) , Re(z) > , Re(fc) > . (A4) 

/o 



t 



with both Re(z) and Re(fc) — »■ 0. With these regulariza- 
tions, the boundary term vanishes and one finds finally 



«A,fe, = 



?rA 



[47rcjfc smli^A)] 1 / 2 



-iX 



/3aV 



= -er^a R K 



(A5) 



Comparing these expressions with (47) yields (50). 

From the definition (Al) and the mode expansions (2) 
and (14) one obtains the relation between the Minkowski 
and Rindler creation and annihilators from which eq. 
(30) follows (one could also obtain them from (26) and 
eq. (A6) below) 

The coefficients A and B in (18) are then easely cal- 
culated. Let us derive the coefficients i^. Substitute 
the expansion (Al) in the r.h.s. of eq. (17a) and regroup 
the terms multiplying (p M and ip M * 



f +00 

ip% = I dk z (a n aS kM +PnPS-k t )<Pk t 



+ (aa8&.+0nag* fc >£* (A6) 
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We used the property (A2). The terms <p^f * are absent 
from (18). Since the Minkowski modes form an orthonor- 
mal basis, each term in the brackets multiplying <p^f * 
must vanish and therefore 



0? 



ilk- 



R* 



a 



(A7) 



Finally substituting these expressions into the first line 
of (A6) one gets with the help of the unitarity relation 
3o| 2 - 1 



\an\ 2 - 



an 



(A8) 



Similarly one obtains = Pnk I 

It is a good check to proove with these expressions the 
orthogonality and completeness of the Unruh modes. For 
instance with and Q! > 0, 

/oo 
dk z A* nK A Q >k z 
-OO 

e w(n+n')/2 



-oo 

1 



anUiv 4^ (sinhTrrisinhTrri') 



1/2 



+oc dh 

-oo w fc 



- K 



j(n-n') 



(A9) 



In the first line we used the orthogonality of Minkowski 
modes (3). The integral becomes trivial after the change 

of variable x = In ) and the factors combines to 

give 

(y>K, "0 ■ (A10) 



Appendix B: Taylor expansions 

We consider dispersion relations admitting a Taylor 
expansion 



n=0 v 



(Bl) 



We assume its convergence radius infinite. We renamed 
«o = i> to be in keeping with the notation (59). We 
substitute (Bl) into the integral expression (69) of the 
Wightman function and expand in terms of the a n >i, 
keeping the first term v in the phase 



y/ — l _ f°° dk ( e -ik(vt+r) c -ik(vt-r)^ c -iktf(k) c -ek 

8n 2 r Jo v + f \ 



J- J^k (e~^ - e~^) e~* f>(*) (£)" . (B2) 



Provided that v > 1, which we will assume in the follow- 
ing, we can replace e~ ek by t — ie, see the discussion of 
sec. II B. Note that even if the DR is polynomial, the 
Taylor expansion of the two-point function has an infinite 
number of terms. We want to integrate this series term- 
by-term over k. The exchange of the series and integral 
is permitted if and only if the series has an infinite con- 
vergence radius. Otherwise the result of the integration 
gives an asymptotic expansion of the Wightman function. 
We were not able to determine the convergence radius of 
the series under the integral sign, but the form of the 



result indicates that we are doing in fact an asymptotic 
expansion. Carrying out this integration we obtain in- 
deed 

1 1 f Wi W2 1 

W(x u x 2 ) = -__j W0+J _-_ + ...j(B3) 

We note s 2 = — {x\ — x 2 ) 2 the invariant distance, and 
we give the expression of the coefficients in terms of the 
Lorentz factor T — (x\ — x^/s 



w = 



Wi — 



8T 3 v 



v(i + r 2 (v 2 - i)) ' (i + r 2 (« 2 -i)) 3 

w 2 = { a2 «[(r 2 -i) 3 -ii w 2 r 2 (r 2 -i) 2 -5« 4 r 4 (r 2 -i) + i5 V 6 r 6 

-2a\ [(r 2 - i) 3 - 5« 2 r 2 (r 2 - 1) 2 + 55 V 4 r 4 (r 2 - 1) + 45« 6 r 6 ] } 



v 3 (i + r 2 (v 2 -i)Y 



(B4) 



(B5) 
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To simplify the discussion we take v = 1 from now on. On inertial trajectories with Lorentz factor 7, the r.h.s. of 
(B3) is given by 

VMr) = ~ 2 {l + + j^y 2 K7 4 (1 + 0(7- 2 )) - 8«?7 6 (1 + 0(7- 2 ))] + •••} , (B6) 



and r stands for r— ie. The series on the r.h.s. of (B6) has 
an essential singularity at r = 0, so it cannot be equal to 
the Wightman function which is a tempered distribution. 
We thus conjecture that for generic dispersion relations, 
the integration term-by-term of (B2) gives an asymptotic 
expansion of the Wightman function. The expansions 
(B3) and (B6) then make sense only for Ms ^> 1 and 
they must be truncated. (The order of truncation chosen 
to minimize the error depends on the detailled behaviour 
of the w n , which we do not know.) 

Now, we want to substitute this expansion into (53) 
and integrate term-by-term in order to get an expansion 



in powers of E/M . To be consistent with Mr 3> 1, we 
are limited to energies 



E < M . (B7) 



Luckily this is compatible with the condition Et 1 
(the Golden rule limit), necessary for the notion of a 
transition rate to be meaningful, and for the expression 
(53) to be valid. The result is 



R 



In 




4a 



7 3 £ 



Acto 



fE_ 



~f 3 E 



(B8) 



We conclude that provided that all the terms retained 
in the sum are small, the transition rate of an inertial 
observer is insensitive to the deviations from Lorentz in- 
variance. From the first two orders we get the restric- 

'V 3 F ■y 4 F 2 

tions 8ai-j^- <C 1 and 24a 2 3 jp- <C 1, and provided all 
the a n are of the same order, the first condition is also 



valid for all the odd order and the second for all the even 
orders. This means that the boost factor is limited to 
7 < {MjEfl 2 or 7 « (M/E) 1 / 3 (depending on the 
presence or not of odd terms). Either way, this is not 
constraining in practice. 



Matters are different for UA trajectories because w n are time dependent. Let us take a.\ = to simplify the 
calculations. Since the Wightman function is not stationary, we remind that we calculate the mean rate (89). We are 
now integrating 



R±(r) 



_2_ r+00 

J —OO 



E 




^ a 2 a 2 24 cosh 4 (ar — x) — 8 cosh 2 (ar — x) — 1 
sinh 2 {x - ie) V 4M 2 sinh 2 (x - ie) 



M 
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\M) 



(i + o(r- 2 )) 



(B9) 



r 



where the Lorentz factor is F = cosh(ar). The correction 
are small for times 



1 o , . max(£^, a) . 

^ ( r 2 (r) [ M > ) 1 



(BIO) 



Since E/M and a/M cannot be larger than 10 100 in prac- 
tice, this means 



ar < 100 , 



(Bll) 



which is still compatible with the condition 
ar > 2tt , 



(B12) 



for the detector to thermalize with a thermal bath at 
temperature a/27r (the typical circular frequency of the 
particles of the bath) . We can conlude that the transition 
rates R± are thermal at best within the interval defined 
by (Bll) and (B12), which represents a few thermal peri- 
odes only (but increasing E or a lowers the upper bound) . 
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Before this, the detector has not thermalized, and beyond 
it, stationarity is lost. 

Our final comment is on [19] in which the transition 
rates were calculated with the help of such a Taylor ex- 
pansion. The author claims to have proved the robust- 
ness of the Unruh effect. He did not. We agree with 
his results, eqs. (25) and (31) (in respectively two and 
four dimensions), obtained from the expression (53) as 
a starting point. As the author noticed, his eqs. (25) 
and (31) are valid only for times r < a -1 . The author 
failed however to recognize that the expression (53) is 
valid precisely in the opposite limit ar > 1. His calcu- 
lations therefore do not prove anything about the Unruh 
effect. 



Appendix C: Steepest descent approximation of eq. 

(102) 

The saddle points solutions of f'(x±) — are given in 
terms of x± — ln(y±) with 



V± 



-r]±VA 
uj — k z 
A = l + k 2 -uj 2 . 



(CI) 



At these points the function / and its derivatives take 
the values 



f(x±) 



i (r]x± ± V~Aj , 
±n/A, f"'(x ± ) 



-IT] . 



(C2) 



The saddle point approximation is good provided |A| » 
1, that is ui » k z and co <C k z . 

We begin with the description of the positions of the 
saddle points in the plane (uj,k z ). The latter is divided 
into three regions: I for < uj < k Zl II for k 2 < uj 2 < 
k 2 + 1, and III for lu 2 > k 2 + 1. Keeping k z fixed and 
increasing uj from to infinity, the saddle points migrate 
as follows. 

In I, A > and y ± ^ 0, hence 



x± 



In 



At?? 

k z — uj 



(C3) 



The real part of (x + , rj - 
decreases from +oo to In 



— 1) and (x_,?7 
y/k* + 1 + fc z ) /k 2 



-1) 
and 



the real part the other two increases from In k z to 



hi 



1 + k z J jk z 

Finally in III, A < and the saddle points are complex 



x± = In 



uj + k z 



=F T)i arctan v 7 — A 



(C6) 



For each process r], both saddle points start at ln(fc 2 + 
i/l + k"l), one migrates to itt/2 and the other to —in/2. 

We now wish to know which of the saddle points are 
actif for a given value of (to,k z ). The steepest descent 
path 7(A), which passes through the actif saddle point(s), 
is by definition a path of stationary phase, i.e. 



Im(/[ 7 (A)])=Im(/(x ± )), 



(C7) 



(or equivalently such that Rc(/(x)) is maximum on 7), 
such that the rate at which Re(/(x)) decreases away from 
the saddle point is as high as possible, that is in a neigh- 
bourhood of the saddle point one requires 



^argtnxiXx^-Xi) 2 ) 



(C8) 



The steepest descent path is a solution of both (C7) and 
(C8). Letting f"(x±) = \f±\e i2a and taking the origin of 
the parameter A at the saddle point, x — x± = \e 1Vsdp + 
0(X 2 ), the angle of the tangent of 7 is given by 



a mod(7r) , 



(C9) 



at the saddle point. In regions I and II, a = ±j (for x± 
respectively, see eq. (C2)), hence if the steepest descent 
path exists, its tangent at the saddle point is (p s dp = ±f ■ 
In region III, a = \ and tp s( i p = 0. 

We give the curves solutions of (C7) in the parametric 
form 7(A) = A + iy(X) and the corresponding angles $ 



/) 7 = A + i arccos 



A 
A 



Rc(x±) 
Re(x±) + 



The real part of (x + ,rf = -1) and (x-,rf = +1) increases II) 7 = A + iarccos 



from In 



(ys/k 2 + 1 + 1^ /k z to +00, and the real part of 



^In^gJ^ix/A-r/A' 
uj sinh A — k z cosh A 

i?± = ±tt ; 
i?± =0 

In ii 2 ^ ± 7A - V X 
k z cosh A 



oj — k z 

uj sinh A 



(v^FTI + i) /k z 



the other two increases from — In 
In k z . 

In II, A > and y± is of the sign of —rj, so we have 



to 



, ,itVA, . 

T/ = +1 , X± = 111 I ; I + ITT 



x± 




III) 



(C4) 
(C5) 



A -4- Rc(.t ± ) , z? ± 
A Rc(.t±)+ , d± 
( 

7 = A ± i arccos 

L 

\ 

i?± =tt(1±1)/2 



tt(»7±1)/2 ; 
n(±l- V )/2 

'/(Iha/W-A 



sinh A — k z cosh A 



(CIO) 

(Cll) 
(C12) 

(C13) 
(C14) 
(C15) 



Along the curves of constant phase (CIO) and (C12) in 
rcsp. regions I and II, 1? changes by ir as the saddle point 
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is passed, so the steepest descent path does not exist. 
Fig. 6 shows Re (/(#)) in the complex x-plane. Indeed, 
the angle of the tangent of the curves is not continuous 
at the saddle point, so it does not verify (C9). In region 



III, the steepest descent path exists. It is the curve of 
stationary phase (C14) passing through the saddle point 
x+. 
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